Genomic and Phylogenetic Characterisation of SARS-CoV-2 Genomes Isolated in Patients from Lambayeque Region, Peru.

OBJECTIVE: this study aims to identify and characterise genomic and phylogenetically isolated SARS-CoV-2 viral isolates in patients from Lambayeque, Peru. METHODS: Nasopharyngeal swabs were taken from patients from the Almanzor Aguinaga Asenjo Hospital, Chiclayo, Lambayeque, Peru, which had been considered mild, moderate, and severe cases of COVID-19. Patients had to have tested positive for COVID-19, using a positive RT-PCR for SARS-CoV-2. Subsequently, the SARS-CoV-2 complete viral genome sequencing was carried out using Illumina MiSeq®. The sequences obtained from the sequence were analysed in Nextclade V1.10.0 to assign the corresponding clades, identify mutations in the SARS-CoV-2 genes and perform quality control of the sequences obtained. All sequences were aligned using MAFFT v7.471. The SARS-CoV-2 isolate Wuhan NC 045512.2 was used as a reference sequence to analyse mutations at the amino acid level. The construction of the phylogenetic tree model was achieved with IQ-TREE v1.6.12. RESULTS: It was determined that during the period from December 2020 to January 2021, the lineages s C.14, C.33, B.1.1.485, B.1.1, B.1.1.1, and B.1.111 circulated, with lineage C.14 being the most predominant at 76.7% (n = 23/30). These lineages were classified in clade 20D mainly and also within clades 20B and 20A. On the contrary, the variants found in the second batch of samples of the period from September to October 2021 were Delta (72.7%), Gamma (13.6%), Mu (4.6%), and Lambda (9.1%), distributed between clades 20J, 21G, 21H, 21J, and 21I. CONCLUSIONS: This study reveals updated information on the viral genomics of SARS-CoV-2 in the Lambayeque region, Peru, which is crucial to understanding the origins and dispersion of the virus and provides information on viral pathogenicity, transmission and epidemiology.

Sphingolipid-Based Synergistic Interactions to Enhance Chemosensitivity in Lung Cancer Cells.

Tumor heterogeneity leads to drug resistance in cancer treatment with the crucial role of sphingolipids in cell fate and stress signaling. We analyzed sphingolipid metabolism and autophagic flux to study chemotherapeutic interactions on the A549 lung cancer model. Loaded cells with fluorescent sphingomyelin analog (BODIPY) and mCherry-EGFP-LC3B were used to track autophagic flux and assess cytotoxicity when cells are exposed to chemotherapy (epirubicin, cisplatin, and paclitaxel) together with sphingolipid pathway inhibitors and autophagy modulators. Our cell model approach employed fluorescent sphingolipid biosensors and a Gaussian Mixture Model of cell heterogeneity profiles to map the influence of chemotherapy on the sphingolipid pathway and infer potential synergistic interactions. Results showed significant synergy, especially when combining epirubicin with autophagy inducers (rapamycin and Torin), reducing cell viability. Cisplatin also synergized with a ceramidase inhibitor. However, paclitaxel often led to antagonistic effects. Our mapping model suggests that combining chemotherapies with autophagy inducers increases vesicle formation, possibly linked to ceramide accumulation, triggering cell death. However, the in silico model proposed ceramide accumulation in autophagosomes, and kinetic analysis provided evidence of sphingolipid colocalization in autophagosomes. Further research is needed to identify specific sphingolipids accumulating in autophagosomes. These findings offer insights into potential strategies for overcoming chemotherapy resistance by targeting the sphingolipid pathway.

Clostridium perfringens phospholipase C, an archetypal bacterial virulence factor, induces the formation of extracellular traps by human neutrophils.

Neutrophil extracellular traps (NETs) are networks of DNA and various microbicidal proteins released to kill invading microorganisms and prevent their dissemination. However, a NETs excess is detrimental to the host and involved in the pathogenesis of various inflammatory and immunothrombotic diseases. Clostridium perfringens is a widely distributed pathogen associated with several animal and human diseases, that produces many exotoxins, including the phospholipase C (CpPLC), the main virulence factor in gas gangrene. During this disease, CpPLC generates the formation of neutrophil/platelet aggregates within the vasculature, favoring an anaerobic environment for C. perfringens growth. This work demonstrates that CpPLC induces NETosis in human neutrophils. Antibodies against CpPLC completely abrogate the NETosis-inducing activity of recombinant CpPLC and C. perfringens secretome. CpPLC induces suicidal NETosis through a mechanism that requires calcium release from inositol trisphosphate receptor (IP3) sensitive stores, activation of protein kinase C (PKC), and the mitogen-activated protein kinase/extracellular signal-regulated kinase (MEK/ERK) pathways, as well as the production of reactive oxygen species (ROS) by the metabolism of arachidonic acid. Proteomic analysis of the C. perfringens secretome identified 40 proteins, including a DNAse and two 5´-nucleotidases homologous to virulence factors that could be relevant in evading NETs. We suggested that in gas gangrene this pathogen benefits from having access to the metabolic resources of the tissue injured by a dysregulated intravascular NETosis and then escapes and spreads to deeper tissues. Understanding the role of NETs in gas gangrene could help develop novel therapeutic strategies to reduce mortality, improve muscle regeneration, and prevent deleterious patient outcomes.

Genomic Characterization of Antibiotic-Resistant Campylobacterales Isolated from Chilean Poultry Meat.

Due to the lack of knowledge about Campylobacterales in the Chilean poultry industry, the objective of this research was to know the prevalence, resistance, and genotypes of Campylobacter, Arcobacter and Helicobacter in 382 samples of chicken meat purchased in Valdivia, Chile. The samples were analyzed using three isolation protocols. Resistance to four antibiotics was evaluated by phenotypic methods. Genomic analyses were performed on selected resistant strains to detect resistance determinants and their genotypes. A total of 59.2% of the samples were positive. Arcobacter butzleri (37.4%) was the most prevalent species, followed by Campylobacter jejuni (19.6%), C. coli (11.3%), A. cryaerophilus (3.7%) and A. skirrowii (1.3%). Helicobacter pullorum (14%) was detected by PCR in a subset of samples. Campylobacter jejuni was resistant to ciprofloxacin (37.3%) and tetracycline (20%), while C. coli and A. butzleri were resistant to ciprofloxacin (55.8% and 2.8%), erythromycin (16.3% and 0.7%) and tetracycline (4.7% and 2.8%), respectively. Molecular determinants were consistent with phenotypic resistance. The genotypes of C. jejuni (CC-21, CC-48, CC-49, CC-257, CC-353, CC-443, CC-446 and CC-658) and C. coli (CC-828) coincided with genotypes of Chilean clinical strains. These findings suggest that besides C. jejuni and C. coli, chicken meat could play a role in the transmission of other pathogenic and antibiotic-resistant Campylobacterales.

In Silico Pipeline to Identify Tumor-Specific Antigens for Cancer Immunotherapy Using Exome Sequencing Data.

Tumor-specific antigens or neoantigens are peptides that are expressed only in cancer cells and not in healthy cells. Some of these molecules can induce an immune response, and therefore, their use in immunotherapeutic strategies based on cancer vaccines has been extensively explored. Studies based on these approaches have been triggered by the current high-throughput DNA sequencing technologies. However, there is no universal nor straightforward bioinformatic protocol to discover neoantigens using DNA sequencing data. Thus, we propose a bioinformatic protocol to detect tumor-specific antigens associated with single nucleotide variants (SNVs) or "mutations" in tumoral tissues. For this purpose, we used publicly available data to build our model, including exome sequencing data from colorectal cancer and healthy cells obtained from a single case, as well as frequent human leukocyte antigen (HLA) class I alleles in a specific population. HLA data from Costa Rican Central Valley population was selected as an example. The strategy included three main steps: (1) pre-processing of sequencing data; (2) variant calling analysis to detect tumor-specific SNVs in comparison with healthy tissue; and (3) prediction and characterization of peptides (protein fragments, the tumor-specific antigens) derived from the variants, in the context of their affinity with frequent alleles of the selected population. In our model data, we found 28 non-silent SNVs, present in 17 genes in chromosome one. The protocol yielded 23 strong binders peptides derived from the SNVs for frequent HLA class I alleles for the Costa Rican population. Although the analyses were performed as an example to implement the pipeline, to our knowledge, this is the first study of an in silico cancer vaccine using DNA sequencing data in the context of the HLA alleles. It is concluded that the standardized protocol was not only able to identify neoantigens in a specific but also provides a complete pipeline for the eventual design of cancer vaccines using the best bioinformatic practices. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-022-00084-9.

De novo assembly of a chromosome-level reference genome of the ornamental butterfly Sericinus montelus based on nanopore sequencing and Hi-C analysis.

Sericinus montelus (Lepidoptera, Papilionidae, Parnassiinae) is a high-value ornamental swallowtail butterfly species widely distributed in Northern and Central China, Japan, Korea, and Russia. The larval stage of this species feeds exclusively on Aristolochia plants. The Aristolochia species is well known for its high levels of aristolochic acids (AAs), which have been found to be carcinogenic for numerous animals. The swallowtail butterfly is among the few that can feed on these toxic host plants. However, the genetic adaptation of S. montelus to confer new abilities for AA tolerance has not yet been well explored, largely due to the limited genomic resources of this species. This study aimed to present a chromosome-level reference genome for S. montelus using the Oxford Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C technology. The final assembly was composed of 581.44 Mb with an expected genome size of 619.27 Mb. Further, 99.98% of the bases could be anchored onto 30 chromosomes. The N50 of contigs and scaffolds was 5.74 and 19.12 Mb, respectively. Approximately 48.86% of the assembled genome was suggested to be repeat elements, and 13,720 protein-coding genes were predicted in the current assembly. The phylogenetic analysis indicated that S. montelus diverged from the common ancestor of swallowtails about 58.57-80.46 million years ago. Compared with related species, S. montelus showed a significant expansion of P450 gene family members, and positive selections on eloa, heatr1, and aph1a resulted in the AA tolerance for S. montelus larva. The de novo assembly of a high-quality reference genome for S. montelus provided a fundamental genomic tool for future research on evolution, genome genetics, and toxicology of the swallowtail butterflies.

Overview of the SARS-CoV-2 genotypes circulating in Latin America during 2021.

Latin America is one of the regions in which the COVID-19 pandemic has a stronger impact, with more than 72 million reported infections and 1.6 million deaths until June 2022. Since this region is ecologically diverse and is affected by enormous social inequalities, efforts to identify genomic patterns of the circulating SARS-CoV-2 genotypes are necessary for the suitable management of the pandemic. To contribute to the genomic surveillance of the SARS-CoV-2 in Latin America, we extended the number of SARS-CoV-2 genomes available from the region by sequencing and analyzing the viral genome from COVID-19 patients from seven countries (Argentina, Brazil, Costa Rica, Colombia, Mexico, Bolivia, and Peru). Subsequently, we analyzed the genomes circulating mainly during 2021 including records from GISAID database from Latin America. A total of 1,534 genome sequences were generated from seven countries, demonstrating the laboratory and bioinformatics capabilities for genomic surveillance of pathogens that have been developed locally. For Latin America, patterns regarding several variants associated with multiple re-introductions, a relatively low percentage of sequenced samples, as well as an increment in the mutation frequency since the beginning of the pandemic, are in line with worldwide data. Besides, some variants of concern (VOC) and variants of interest (VOI) such as Gamma, Mu and Lambda, and at least 83 other lineages have predominated locally with a country-specific enrichments. This work has contributed to the understanding of the dynamics of the pandemic in Latin America as part of the local and international efforts to achieve timely genomic surveillance of SARS-CoV-2.

Difference in mortality rates in hospitalized COVID-19 patients identified by cytokine profile clustering using a machine learning approach: An outcome prediction alternative.

COVID-19 is a disease caused by the novel Coronavirus SARS-CoV-2 causing an acute respiratory disease that can eventually lead to severe acute respiratory syndrome (SARS). An exacerbated inflammatory response is characteristic of SARS-CoV-2 infection, which leads to a cytokine release syndrome also known as cytokine storm associated with the severity of the disease. Considering the importance of this event in the immunopathology of COVID-19, this study analyses cytokine levels of hospitalized patients to identify cytokine profiles associated with severity and mortality. Using a machine learning approach, 3 clusters of COVID-19 hospitalized patients were created based on their cytokine profile. Significant differences in the mortality rate were found among the clusters, associated to different CXCL10/IL-38 ratio. The balance of a CXCL10 induced inflammation with an appropriate immune regulation mediated by the anti-inflammatory cytokine IL-38 appears to generate the adequate immune context to overrule SARS-CoV-2 infection without creating a harmful inflammatory reaction. This study supports the concept that analyzing a single cytokine is insufficient to determine the outcome of a complex disease such as COVID-19, and different strategies incorporating bioinformatic analyses considering a broader immune profile represent a more robust alternative to predict the outcome of hospitalized patients with SARS-CoV-2 infection.

Clinical Profiles at the Time of Diagnosis of SARS-CoV-2 Infection in Costa Rica During the Pre-vaccination Period Using a Machine Learning Approach.

The clinical manifestations of COVID-19, caused by the SARS-CoV-2, define a large spectrum of symptoms that are mainly dependent on the human host conditions. In Costa Rica, more than 169,000 cases and 2185 deaths were reported during the year 2020, the pre-vaccination period. To describe the clinical presentations at the time of diagnosis of SARS-CoV-2 infection in Costa Rica during the pre-vaccination period, we implemented a symptom-based clustering using machine learning to identify clusters or clinical profiles at the population level among 18,974 records of positive cases. Profiles were compared based on symptoms, risk factors, viral load, and genomic features of the SARS-CoV-2 sequence. A total of 18 symptoms at time of diagnosis of SARS-CoV-2 infection were reported with a frequency > 1%, and those were used to identify seven clinical profiles with a specific composition of clinical manifestations. In the comparison between clusters, a lower viral load was found for the asymptomatic group, while the risk factors and the SARS-CoV-2 genomic features were distributed among all the clusters. No other distribution patterns were found for age, sex, vital status, and hospitalization. In conclusion, during the pre-vaccination time in Costa Rica, the symptoms at the time of diagnosis of SARS-CoV-2 infection were described in clinical profiles. The host co-morbidities and the SARS-CoV-2 genotypes are not specific of a particular profile, rather they are present in all the groups, including asymptomatic cases. In addition, this information can be used for decision-making by the local healthcare institutions (first point of contact with health professionals, case definition, or infrastructure). In further analyses, these results will be compared against the profiles of cases during the vaccination period. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-022-00058-x.

Metagenomic pipeline for identifying co-infections among distinct SARS-CoV-2 variants of concern: study cases from Alpha to Omicron.

Concomitant infection or co-infection with distinct SARS-CoV-2 genotypes has been reported as part of the epidemiological surveillance of the COVID-19 pandemic. In the context of the spread of more transmissible variants during 2021, co-infections are not only important due to the possible changes in the clinical outcome, but also the chance to generate new genotypes by recombination. However, a few approaches have developed bioinformatic pipelines to identify co-infections. Here we present a metagenomic pipeline based on the inference of multiple fragments similar to amplicon sequence variant (ASV-like) from sequencing data and a custom SARS-CoV-2 database to identify the concomitant presence of divergent SARS-CoV-2 genomes, i.e., variants of concern (VOCs). This approach was compared to another strategy based on whole-genome (metagenome) assembly. Using single or pairs of sequencing data of COVID-19 cases with distinct SARS-CoV-2 VOCs, each approach was used to predict the VOC classes (Alpha, Beta, Gamma, Delta, Omicron or non-VOC and their combinations). The performance of each pipeline was assessed using the ground-truth or expected VOC classes. Subsequently, the ASV-like pipeline was used to analyze 1021 cases of COVID-19 from Costa Rica to investigate the possible occurrence of co-infections. After the implementation of the two approaches, an accuracy of 96.2% was revealed for the ASV-like inference approach, which contrasts with the misclassification found (accuracy 46.2%) for the whole-genome assembly strategy. The custom SARS-CoV-2 database used for the ASV-like analysis can be updated according to the appearance of new VOCs to track co-infections with eventual new genotypes. In addition, the application of the ASV-like approach to all the 1021 sequenced samples from Costa Rica in the period October 12th-December 21th 2021 found that none corresponded to co-infections with VOCs. In conclusion, we developed a metagenomic pipeline based on ASV-like inference for the identification of co-infection with distinct SARS-CoV-2 VOCs, in which an outstanding accuracy was achieved. Due to the epidemiological, clinical, and molecular relevance of the concomitant infection with distinct genotypes, this work represents another piece in the process of the surveillance of the COVID-19 pandemic in Costa Rica and worldwide.

Analysis of gut microbiome profiles in common marmosets (Callithrix jacchus) in health and intestinal disease.

Chronic gastrointestinal (GI) diseases are the most common diseases in captive common marmosets. To understand the role of the microbiome in GI diseases, we characterized the gut microbiome of 91 healthy marmosets (303 samples) and 59 marmosets diagnosed with inflammatory bowel disease (IBD) (200 samples). Healthy marmosets exhibited "humanized," Bacteroidetes-dominant microbiomes. After up to 2 years of standardized diet, housing and husbandry, marmoset microbiomes could be classified into four distinct marmoset sources based on Prevotella and Bacteroides levels. Using a random forest (RF) model, marmosets were classified by source with an accuracy of 93% with 100% sensitivity and 95% specificity using abundance data from 4 Prevotellaceae amplicon sequence variants (ASVs), as well as single ASVs from Coprobacter, Parabacteroides, Paraprevotella, Phascolarctobacterium, Oribacterium and Fusobacterium. A single dysbiotic IBD state was not found across all marmoset sources, but IBD was associated with lower alpha diversity and a lower Bacteroides:Prevotella copri ratio within each source. IBD was highest in a Prevotella-dominant cohort, and consistent with Prevotella-linked diseases, pro-inflammatory genes in the jejunum were upregulated. RF analysis of serum biomarkers identified serum calcium, hemoglobin and red blood cell (RBC) counts as potential biomarkers for marmoset IBD. This study characterizes the microbiome of healthy captive common marmosets and demonstrates that source-specific microbiomes can be retained despite standardized diets and husbandry practices. Marmosets with IBD had decreased alpha diversity and a shift in the ratio of Bacteroides:Prevotella copri compared to healthy marmosets.

Alterations in common marmoset gut microbiome associated with duodenal strictures.

Chronic gastrointestinal (GI) diseases are the most common diseases in captive common marmosets (Callithrix jacchus). Despite standardized housing, diet and husbandry, a recently described gastrointestinal syndrome characterized by duodenal ulcers and strictures was observed in a subset of marmosets sourced from the New England Primate Research Center. As changes in the gut microbiome have been associated with GI diseases, the gut microbiome of 52 healthy, non-stricture marmosets (153 samples) were compared to the gut microbiome of 21 captive marmosets diagnosed with a duodenal ulcer/stricture (57 samples). No significant changes were observed using alpha diversity metrics, and while the community structure was significantly different when comparing beta diversity between healthy and stricture cases, the results were inconclusive due to differences observed in the dispersion of both datasets. Differences in the abundance of individual taxa using ANCOM, as stricture-associated dysbiosis was characterized by Anaerobiospirillum loss and Clostridium perfringens increases. To identify microbial and serum biomarkers that could help classify stricture cases, we developed models using machine learning algorithms (random forest, classification and regression trees, support vector machines and k-nearest neighbors) to classify microbiome, serum chemistry or complete blood count (CBC) data. Random forest (RF) models were the most accurate models and correctly classified strictures using either 9 ASVs (amplicon sequence variants), 4 serum chemistry tests or 6 CBC tests. Based on the RF model and ANCOM results, C. perfringens was identified as a potential causative agent associated with the development of strictures. Clostridium perfringens was also isolated by microbiological culture in 4 of 9 duodenum samples from marmosets with histologically confirmed strictures. Due to the enrichment of C. perfringens in situ, we analyzed frozen duodenal tissues using both 16S microbiome profiling and RNAseq. Microbiome analysis of the duodenal tissues of 29 marmosets from the MIT colony confirmed an increased abundance of Clostridium in stricture cases. Comparison of the duodenal gene expression from stricture and non-stricture marmosets found enrichment of genes associated with intestinal absorption, and lipid metabolism, localization, and transport in stricture cases. Using machine learning, we identified increased abundance of C. perfringens, as a potential causative agent of GI disease and intestinal strictures in marmosets.

Insights into the mutation T1117I in the spike and the lineage B.1.1.389 of SARS-CoV-2 circulating in Costa Rica.

Emerging mutations and genotypes of the SARS-CoV-2 virus, responsible for the COVID-19 pandemic, have been reported globally. In Costa Rica during the year 2020, a predominant genotype carrying the mutation T1117I in the spike (S:T1117I) was previously identified. To investigate the possible effects of this mutation on the function of the spike, i.e. the biology of the virus, different bioinformatic pipelines based on phylogeny, natural selection, and co-evolutionary models, molecular docking, and epitopes prediction were implemented. Results of the phylogeny of sequences carrying the S:T1117I worldwide showed a polyphyletic group, with the emergence of local lineages. In Costa Rica, the mutation is found in the lineage B.1.1.389 and it is suggested to be a product of positive/adaptive selection. Different changes in the function of the spike protein and more stable interaction with a ligand (nelfinavir drug) were found. Only one epitope out 742 in the spike was affected by the mutation, with some different properties, but suggesting scarce changes in the immune response and no influence on the vaccine effectiveness. Jointly, these results suggest a partial benefit of the mutation for the spread of the virus with this genotype during the year 2020 in Costa Rica, although possibly not strong enough with the introduction of new lineages during early 2021 which became predominant later. In addition, the bioinformatic analyses used here can be applied as an in silico strategy to eventually study other mutations of interest for the SARS-CoV-2 virus and other pathogens.

Avances de la bioinformática en Costa Rica: vista retrospectiva y perspectivas / Bioinformatics advances in Costa Rica: retrospective view and perspectives

Resumen Introducción: La disciplina científica de la bioinformática tiene el potencial de generar aplicaciones innovadoras para las sociedades humanas. Costa Rica, pequeña en tamaño y población en comparación con otros países de América Latina, ha ido adoptando la disciplina de manera progresiva. El reconocer los avances permite determinar hacia dónde puede dirigirse el país en este campo, así como su contribución a la región latinoamericana. Objetivo: En este manuscrito se reporta evidencia de la evolución de la bioinformática en Costa Rica, para identificar debilidades y fortalezas que permitan definir acciones a futuro. Métodos: Se realizaron búsquedas en bases de datos de publicaciones científicas y repositorios de secuencias, así como información de actividades de capacitación, redes, infraestructura, páginas web y fuentes de financiamiento. Resultados: Se observan avances importantes desde el 2010, incluyendo un aumento en oportunidades de entrenamiento y número de publicaciones, aportes significativos a las bases de datos de secuencias y conexiones por medio de redes. Sin embargo, ciertas áreas, como la masa crítica y la financiación requieren más desarrollo. La comunidad científica y sus patrocinadores deben promover la investigación basada en bioinformática, invertir en la formación de estudiantes de posgrado, aumentar la formación de profesionales, crear oportunidades laborales para carreras en bioinformática y promover colaboraciones internacionales a través de redes. Conclusiones: Se sugiere que para experimentar los beneficios de las aplicaciones de la bioinformática se deben fortalecer tres aspectos clave: la comunidad científica, la infraestructura de investigación y las oportunidades de financiamiento. El impacto de tal inversión sería el desarrollo de proyectos ambiciosos pero factibles y colaboraciones extendidas dentro de la región latinoamericana. Esto permitiría realizar contribuciones significativas para abordar los desafíos globales y la aplicación de nuevos enfoques de investigación, innovación y transferencia de conocimiento para el desarrollo de la economía, dentro de un marco de ética de la investigación.

Abstract Introduction: The scientific discipline of bioinformatics has the potential to generate innovative applications for human societies. Costa Rica, small in size and population compared to other Latin American countries, has been progressively adopting the discipline. Recognizing progress makes it possible to determine where the country can go in this field, as well as its contribution to the Latin American region. Objective: This manuscript reports evidence of the evolution of bioinformatics in Costa Rica, to identify weaknesses and strengths allowing future actions plans. Methods: We searched databases of scientific publications and sequence repositories, as well as information on training activities, networks, infrastructure, web pages and funding sources. Results: Important advances have been observed since 2010, such as increases in training opportunities and the number of publications, significant contributions to the sequence databases and connections through networks. However, areas such as critical mass and financing require further development. The scientific community and its sponsors should promote bioinformatics-based research, invest in graduate student training, increase professional training, create career opportunities in bioinformatics, and promote international collaborations through networks. Conclusions: It is suggested that in order to experience the benefits of bioinformatics applications, three key aspects must be strengthened: the scientific community, the research infrastructure, and funding opportunities. The impact of such investment would be the development of ambitious but feasible projects and extended collaborations within the Latin American region and abroad. This would allow significant contributions to address global challenges and the implementation of new approaches to research, innovation and knowledge transfer for the development of the economy, within an ethics of research framework.

High Efficacy of Therapeutic Equine Hyperimmune Antibodies Against SARS-CoV-2 Variants of Concern.

SARS-CoV-2 variants of concern show reduced neutralization by vaccine-induced and therapeutic monoclonal antibodies; therefore, treatment alternatives are needed. We tested therapeutic equine polyclonal antibodies (pAbs) that are being assessed in clinical trials in Costa Rica against five globally circulating variants of concern: alpha, beta, epsilon, gamma and delta, using plaque reduction neutralization assays. We show that equine pAbs efficiently neutralize the variants of concern, with inhibitory concentrations in the range of 0.146-1.078 µg/mL, which correspond to extremely low concentrations when compared to pAbs doses used in clinical trials. Equine pAbs are an effective, broad coverage, low-cost and a scalable COVID-19 treatment.

SARS-CoV-2 genomic surveillance in Costa Rica: Evidence of a divergent population and an increased detection of a spike T1117I mutation.

Genome sequencing is a key strategy in the surveillance of SARS-CoV-2, the virus responsible for the COVID-19 pandemic. Latin America is the hardest-hit region of the world, accumulating almost 20% of COVID-19 cases worldwide. In Costa Rica, from the first detected case on March 6th to December 31st almost 170,000 cases have been reported. We analyzed the genomic variability during the SARS-CoV-2 pandemic in Costa Rica using 185 sequences, 52 from the first months of the pandemic, and 133 from the current wave. Three GISAID clades (G, GH, and GR) and three PANGOLIN lineages (B.1, B.1.1, and B.1.291) were predominant, suggesting multiple re-introductions from other regions. The whole-genome variant calling analysis identified a total of 283 distinct nucleotide variants, following a power-law distribution with 190 single nucleotide mutations in a single sequence, and only 16 mutations were found in >5% sequences. These mutations were distributed through the whole genome. The prevalence of worldwide-found variant D614G in the Spike (98.9% in Costa Rica), ORF8 L84S (1.1%) is similar to what is found elsewhere. Interestingly, the frequency of mutation T1117I in the Spike has increased during the current pandemic wave beginning in May 2020 in Costa Rica, reaching 29.2% detection in the full genome analyses in November 2020. This variant has been observed in less than 1% of the GISAID reported sequences worldwide in 2020. Structural modeling of the Spike protein with the T1117I mutation suggests a potential effect on the viral oligomerization needed for cell infection, but no differences with other genomes on transmissibility, severity nor vaccine effectiveness are predicted. In conclusion, genome analyses of the SARS-CoV-2 sequences over the course of the COVID-19 pandemic in Costa Rica suggest the introduction of lineages from other countries and the detection of mutations in line with other studies, but pointing out the local increase in the detection of Spike-T1117I variant. The genomic features of this virus need to be monitored and studied in further analyses as part of the surveillance program during the pandemic.

A first perturbome of Pseudomonas aeruginosa: Identification of core genes related to multiple perturbations by a machine learning approach.

Tolerance to stress conditions is vital for organismal survival, including bacteria under specific environmental conditions, antibiotics, and other perturbations. Some studies have described common modulation and shared genes during stress response to different types of disturbances (termed as perturbome), leading to the idea of central control at the molecular level. We implemented a robust machine learning approach to identify and describe genes associated with multiple perturbations or perturbome in a Pseudomonas aeruginosa PAO1 model. Using microarray datasets from the Gene Expression Omnibus (GEO), we evaluated six approaches to rank and select genes: using two methodologies, data single partition (SP method) or multiple partitions (MP method) for training and testing datasets, we evaluated three classification algorithms (SVM Support Vector Machine, KNN K-Nearest neighbor and RF Random Forest). Gene expression patterns and topological features at the systems level were included to describe the perturbome elements. We were able to select and describe 46 core response genes associated with multiple perturbations in P. aeruginosa PAO1 and it can be considered a first report of the P. aeruginosa perturbome. Molecular annotations, patterns in expression levels, and topological features in molecular networks revealed biological functions of biosynthesis, binding, and metabolism, many of them related to DNA damage repair and aerobic respiration in the context of tolerance to stress. We also discuss different issues related to implemented and assessed algorithms, including data partitioning, classification approaches, and metrics. Altogether, this work offers a different and robust framework to select genes using a machine learning approach.

Set-theory based benchmarking of three different variant callers for targeted sequencing.

BACKGROUND: Next generation sequencing (NGS) technologies have improved the study of hereditary diseases. Since the evaluation of bioinformatics pipelines is not straightforward, NGS demands effective strategies to analyze data that is of paramount relevance for decision making under a clinical scenario. According to the benchmarking framework of the Global Alliance for Genomics and Health (GA4GH), we implemented a new simple and user-friendly set-theory based method to assess variant callers using a gold standard variant set and high confidence regions. As model, we used TruSight Cardio kit sequencing data of the reference genome NA12878. This targeted sequencing kit is used to identify variants in key genes related to Inherited Cardiac Conditions (ICCs), a group of cardiovascular diseases with high rates of morbidity and mortality. RESULTS: We implemented and compared three variant calling pipelines (Isaac, Freebayes, and VarScan). Performance metrics using our set-theory approach showed high-resolution pipelines and revealed: (1) a perfect recall of 1.000 for all three pipelines, (2) very high precision values, i.e. 0.987 for Freebayes, 0.928 for VarScan, and 1.000 for Isaac, when compared with the reference material, and (3) a ROC curve analysis with AUC > 0.94 for all cases. Moreover, significant differences were obtained between the three pipelines. In general, results indicate that the three pipelines were able to recognize the expected variants in the gold standard data set. CONCLUSIONS: Our set-theory approach to calculate metrics was able to identify the expected ICCs related variants by the three selected pipelines, but results were completely dependent on the algorithms. We emphasize the importance to assess pipelines using gold standard materials to achieve the most reliable results for clinical application.

Genomic context of the two integrons of ST-111 Pseudomonas aeruginosa AG1: A VIM-2-carrying old-acquaintance and a novel IMP-18-carrying integron.

Pseudomonas aeruginosa is an opportunist and versatile organism responsible for infections mainly in immunocompromised hosts. This pathogen has high intrinsic resistance to most antimicrobials. P. aeruginosa AG1 (PaeAG1) is a Costa Rican high-risk ST-111 strain with resistance to multiple antibiotics, including carbapenems, due to the activity of VIM-2 and IMP-18 metallo-ß-lactamases (MBLs). These genes are harbored in two class 1 integrons located inone out of the 57 PaeAG1 genomic islands. However, the genomic context associated to these determinants in PaeAG1 and other P. aeruginosa strains is unclear. Thus, we first assessed the transcriptional activity of VIM-2 and IMP-18 genes when exposed to imipenem (a carbapenem) by RT-qPCR. To select related genomes to PaeAG1, we implemented a pan-genome analysis to define and up-date the phylogenetic relationship among complete P. aeruginosa genomes. We also studied the PaeAG1 genomic islands content in the related strains and finally we described the architecture and possible evolutionary steps of the genomic regions around the VIM-2- and IMP-18-carrying integrons. Expression of VIM-2 and IMP-18 genes was demonstrated to be induced after imipenem exposure. In a subsequent comparative genomics analysis with 211 strains, the P. aeruginosa pan-genome revealed that complete genome sequences are able to separate clones by MLST profile, including a clear ST-111 cluster with PaeAG1. The PaeAG1 genomic islands were found to define a diverse presence/absence pattern among related genomes. Finally, landscape reconstruction of genomic regions showed that VIM-2-carrying integron (In59-like) is an old-acquaintance element harbored in the same known region found in other two ST-111 strains. Also, PaeAG1 has an exclusive genomic region containing a novel IMP-18-carrying integron (registered as In1666), with an arrangement never reported before. Altogether, we provide new insights about the genomic determinants associated with the resistance to carbapenems in this high-risk P. aeruginosa using comparative genomics.